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Abstract 

We make several observations concerning the low quark mass region with Wilson 
fermions and how this is connected with the e regime in the continuum. A transi- 
tion from tiny cutoff effects to rather large discretization errors would take place in 
general with Wilson fermions if we lower the quark mass at finite lattice spacing. 
We argue that these two regions exhibit rather different behaviours concerning the 
coupling between cutoff effects and zero-modes. We interpolate between these two 
regimes adding to the continuum e regime formulas , in the spirit of the Symanzik ex- 
pansion, the relevant operators parametrising the leading cutoff effects. We compute 
the partition function, the chiral condensate, scalar and pseudo-scalar correlation 
functions. The final formulae can be used to fit lattice data to extract physical low 
energy constants, and to estimate systematic uncertainties coming from discretiza- 
tion errors. Moreover they suggest ways on how to remove these cutoff effects, the 
core of which are captured by the continuum zero modes. 
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1 Introduction 



Simulations of lattice QCD in the so-called e regime [1] of the chiral expansion allow 
in principle the extraction of physical parameters, like decay constants and electroweak 
effective couplings. Up to now almost all the published results of dynamical simulations 
in the e regime have been always obtained using lattice chiral invariant formulations, 
like overlap [2] (and references therein), and perfect action [3] fermions. There are many 
reasons why chiral invariant lattice fermions are preferable. Here we mention the possi- 
bility of lowering the quark mass at finite lattice spacing, without encountering stability 
or metastabilities problems, small cutoff effects for a wide range of quark masses, 0(a 2 ) 
scaling violations and continuum-like renormalization patterns. 

Contrary to what is widely believed simulations in the e regime are in principle not 
exclusive to lattice actions with exact lattice chiral symmetry. It is thus interesting to 
understand how the continuum e regime is probed by Wilson fermions, by mean of an- 
alytical and numerical tools. It has been stressed in the past [4,5] that topology plays 
an important role in this extreme regime. One could still think of a strategy where sim- 
ulations performed with Wilson-like fermions would sample all the topological sectors. 
This could be an alternative way to the e regime. Preliminary and encouraging results of 
dynamical simulations in the e regime using Wilson twisted mass (Wtm) fermions have 
been presented at the last lattice conferences [6,7]. Recently it has been shown [8] that 
a suitable algorithm can sample the configuration space of the e regime also with stan- 
dard Wilson fermions. In order to properly interpret the results of these simulations, it 
is important to understand how Wilson-like fermions probe the low quark mass region at 
finite volume. It is important to understand analytically the quark mass, lattice spacing 
and volume dependence with Wilson fermions. In this Letter we perform a step towards 
this understanding. Our analysis will concern Wilson fermions, but it can be extended to 
Wilson twisted mass. 

It is well known [1] that in continuum QCD if one sends the quark mass to zero keeping 
the size of the volume fixed, chiral symmetry is restored, and the dependence on the quark 
mass of the chiral condensate is smooth. It is also well known that at finite lattice spacing 
and infinite volume Wilson fermions exhibit a peculiar chiral phase diagram [9,10] which 
has important consequences for the quark mass dependence of the condensates and of the 
would-be continuum Goldstone bosons. 

It is interesting and important to understand which mechanism takes place at finite lattice 
spacing, if we lower the quark mass keeping the size of the volume fixed. One way to study 
this is to include the effects of the non- vanishing lattice spacing in the analysis of [1], to 
address these issues from the chiral effective theory point of view. 

The Letter is organized as follows: in sect. 2 we discuss the different patterns of symmetry 
breaking in the continuum and at finite lattice spacing and the connection with the 
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choice of the power counting. This will bring us to find a phenomenological prescription 
on how to include cutoff effects in our computations. Using this prescription in sect. 3 we 
analyze and compute partition function and chiral condensate, and in sect. 4 the scalar 
and pseudoscalar correlation functions. A discussion of the results and conclusions will be 
given in sect. 5. 



2 Power counting 

The e expansion in chiral perturbation theory was introduced in [1] and further developed 
in [11,12,13]. It is interesting to notice that in these papers no mention is made on the 
topological sectors of the configuration space of the theory. In this Letter we will only 
consider the N = 2 flavours case. 

The basic feature of the e expansion is that while the zero modes are treated non- 
perturbatively, the non-zero modes fluctuate in a Gaussian way around them and are 
treated in a standard perturbative way. This is needed to describe correctly the collective 
behaviour of the zero modes when the quark mass is sent to zero at fixed finite volume. 
To obtain this in an algebraic way the standard p power counting is modified treating the 
would be Goldstone boson mass small compared with the box size. The standard power 
counting in the e regime in the continuum reads 



where m indicates the quark mass and L and T the spatial and temporal extent of the 
volume V = L 3 x T. 

It is pedagogical to understand which are the lattice parameters which ought to be used 
to reach this regime. Let us take as a typical range of volumes for dynamical simulations 
with Wilson fermions 1.5fm < L < 3fm. The e regime is reached if mHV < 1 . If we take 
the reasonable value of £ = (250MeV) 3 one easily obtains that, if T = L, the value of the 
properly renormalized quark mass should be of the order of 

m~15 Mev for L = 1.5 fm (2.2) 
m~6 Mev for L = 2 fm (2.3) 
m ~ 2 Mev for L = 3 fm. (2.4) 

While this is just an order of magnitude estimate, it is clear that quark masses significantly 
lower than 20 MeV are needed to be in the e regime with the current available volumes. 

It is interesting now to understand, given these values, and given the current available 
lattice spacings which is correct power counting to adopt in describing the lattice data. In 




(2.1) 
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fact one can include the effect of the discretization errors in the chiral effective theory [10] 
using as a starting point the continuum Symanzik action, describing the interactions of 
quarks and gluons with momenta much smaller than n/a. The symmetry properties of 
the Symanzik action allow to write a generalized chiral expansion with explicit factors of 
the lattice spacing a. This has been done in the past in the p regime with great success for 
a set of lattice actions [14,15,16,17] and for different power countings. The most widely 
used power countings are defined by the way the quark mass and the lattice spacing 
are related to each other. The so-called GSM (generically small masses) regime [17] is 
defined by m ~ a A 2 while the Aoki regime [14], also called large cutoff effects regime 
(LCE) [18], is defined by m ~ a 2 A 3 . For illustrative purposes let us take the reasonable 
value A = 250MeV. Nowadays dynamical lattice simulations with Wilson fermions are 
performed in the range 0.08fm < a < 0.04fm. This gives us the following values 



aA 2 ~ 25 MeV a 2 A 3 ~ 3 MeV for a = 0.08 fm (2.5) 
aA 2 ~ 12 MeV a 2 A 3 ~ 1 MeV for a = 0.04 fm. (2.6) 

A good simulation setup would be with a = 0.04 fm, a lattice box of L/a = T/a = 48 
and a quark mass m = 6 MeV. In this rather ideal case we would be in a region of quark 
masses in the GSM regime or quite close to it. Certainly not everyone will have access 
to such high quality gauge configurations, and moreover it could be advisable anyhow to 
have more lattice spacings to check for scaling violations. If we would have to increase 
the value of the lattice spacing, we would move towards the LCE region possibly getting 
dangerously close to it, In fact in the less ambitious setup where a = 0.08 fm, the lattice 
box is L/a = T/a = 24 and the quark mass m = 6 MeV we would be quite close to 
the LCE region. Even if this discussion is quite simple it gives the hint that in the near 
future if we want to probe the e regime with Wilson fermions, we will simulate in the 
GSM regime but going closer to the LCE regime 

It is thus very important to understand the behaviour of Wilson fermions for small quark 
masses in a finite volume. 



2.1 Generic small masses 



In infinite volume the GSM power counting is such that m and a are of the same order, 
as p 2 . With this power counting the leading order (LO) chiral Lagrangian is 



C 



(2) 

w x 



Tr [9 /1 C/(x) + 9 M C/'(z)J -2B Tr [M*U(x) + MU{x) ] \ -2aW Tr [U(x) + U(x) ] 

(2.7) 

where F and Bq are the LECs appearing at LO, and Wo is an unknown dimensionful 
constant which parametrises the leading cutoff effects, a being the lattice spacing. M. is 
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the quark mass matrix and U is the field collecting the Goldstone bosons fields. We recall 
that if the csw coefficient would be set to its "correct" non-perturbative value we would 
have Wq = 0. This does not mean that all the 0(a) terms will disappear, because there 
will be 0(a) terms at higher orders in the chiral expansion, that would have to be, if 
needed, cancelled by other improvement coefficients. The leading 0(a) can be reabsorbed 
in the definition of the quark mass [10] 

Wn 

M^M' = M + -^a. (2.8) 

-DO 

Thus with this power counting the leading order (LO) Lagrangian is identical to the 
continuum LO Lagrangian. 

Let us briefly recall why one integrates exactly over the zero modes in the e regime in the 
continuum. If one considers the tree-level pion propagator 

G{x) = v^ WTm! = VmI + v J ^TMf (2 ' 9) 

one sees immediately that the zero mode contribution explodes in the chiral limit at fixed 
finite volume. The standard p expansion would fail in this situation because it treats 
all the modes in the same way. Alternatively one can say that in the massless limit the 
zero modes contribution does not appear in the quadratic approximation of the action. To 
circumvent this problem one should reorder the expansion in order to have a representation 
of the chiral effective theory with an arbitrary number of zero modes propagators. This is 
achieved treating the constant zero modes in a collective non-perturbative way, and the 
non-zero modes as standard perturbations [1]. 

This tells us that at finite lattice spacing in the GSM regime the LO behaviour of the zero 
modes is like in the continuum, and to correctly include the zero modes propagation in 
the computation one would have to treat them in the same way, i.e. the power counting 
for GSM in the e regime is 

m = 0(e% | = 0(6), ^ = 0(e) a = 0(e% (2.10) 

To understand the possible modifications induced by the discretization errors at NLO we 
have to scrutinize the NLO chiral Lagrangian of the GSM regime given for example in 
refs. [14,17]. The NLO Lagrangian has a continuum part and a part coming from the 
discretization effects which reads 



d§ = + aWTr^Wd^TriU + U ] ) - 2aB WTi{M' ] U + U ] M')Tr(U + U ] ) 



-w x - '-x 
-a 2 W 



Tr(U + U ] )\ - 2aB H'Tr(M' + M' ] ). (2.11) 
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The notation used here is the same used in ref. [17], but here we redefine the LECs 
absorbing the Wo term. As a result their dimension is different and given by W,W = 

Energy 3 , and [W] = Energy 6 . We also do not include all the terms needed for the 
analysis of the axial and vector currents which we do not consider in this Letter. The 
terms of the kind 0(am) (W) and 0(a 2 ) (W) are of order 0(e 4 )[I]. The same is true for 
the terms of 0(ap 2 ) (W). 

To explain this let us consider the standard parametrisation for the Goldstone field U in 
the e expansion 

r m(x) 



U (x) = U exp 

where ir = n a T a and the r are normalized such that {r a ,r b } = 25 ab . 



(2.12) 



In the e expansion the fluctuations around the constant zero modes are treated as 0(e), 
thus each power of the derivative (or of the momentum p) carries at least one power of 
the fluctuation. As a result the 0(ap 2 ) terms are also of 0(e 4 ). 

The conclusion is that in the GSM regime the e expansion with Wilson fermions is like 
in the continuum up to 0(e 4 ), i.e. if we are in a region of parameters where we can 
neglect NNLO terms, Wilson fermions are automatically 0(a) improved. This is not so 
surprising for the following reason. In the e expansion in the continuum the 0(e°) and 0(e 2 ) 
contributions can be computed considering only the LO Lagrangian. In this particular 
region contributions from what are called NLO LECs' in the p regime, are suppressed by 
one order compared to the standard p regime expansion. 



2.2 Large cutoff effects region 



If we now consider coarser lattice spacings we enter in the LCE regime and the power 
counting changes such that m! and a 2 terms are of the same order, as p 2 . Already at 
leading order the situation is different. The LO Lagrangian is 



C 



(2) 

w x 



— Tr 



d^UixYd^Uix) 



-^Tr 
2 



M n U(x) + M'U(xY] -a 2 W 



Tr 



U(x) + u{ x y 

(2.13) 

where W is a LEC which parametrises 0(a 2 ) cutoff effects. It particular its value depends 
on the value of csw adopted but it does not vanish if the theory is non-perturbatively 
improved. In infinite volume the Lagrangian (2.13) implies a competition of the mass 
term and of the 0(a 2 ) term in the shape of the potential that causes a non-trivial vacuum 
structure. Minimizing the potential gives rise to two possible scenarios [10,16,19,20] for 



1 In the e expansion one has to consider the action including the space-time integration which 
gives a volume factor of 0(e~ 4 ). 
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the phase diagram of Wilson-like fermions: 



• the Aoki scenario W < [9]; 

• the Sharpe-Singleton scenario W > [10]. 

Depending on the sign of W the two scenarios predict a different pattern for the quark 
mass dependence of the chiral condensate and of the pion masses. In the Sharpe-Singleton 
scenario the pion masses do not go to zero in the chiral limit, while in the Aoki scenario 
within the Aoki phase the charged pion stay massless while the neutral pion mass becomes 
massive. Recalling our discussion in the previous section this means that for the Sharpe- 
Singleton scenario the pion mass will never vanish, and as a consequence one will still 
have the zero-modes in the quadratic approximation. This implies that in this scenario 
a standard perturbative expansion can be performed and no collective phenomena take 
place [21]. In the Aoki scenario, the Goldstone boson manifold is not SU{2) like in the 
continuum, and an appropriate expansion should be performed taking into account the 
pattern of symmetry breaking at finite lattice spacing [21]. This does not mean that the 
e regime does not exist for Wilson fermions. In fact for example if we stay in the GSM 
regime, like we have discussed in the previous section, one could simulate in the e regime 
being affected by very tiny cutoff effects. It simply means that the order of the chiral and 
continuum limit has to be understood with care. 

From this discussion it is clear that Wilson fermions in the low quark mass region regime 
have a transition in terms of cutoff effects from no cutoff effects (up to NNLO) in the 
GSM regime to a regime where cutoff effects appear at LO and the collective zero-modes 
phenomena, if present, are pretty different from the continuum. 

To understand up to which value of the lattice spacing and quark mass we can still 
discuss of continuum zero modes we can imagine a situation where we have continuum- 
like collective phenomena and the Lagrangian of the LCE regime (2.13). While we know 
this is not a completely consistent procedure it might give us a hint to where this picture 
breaks down. 

The power counting would be 



and we expect to encounter some contradictions, indicating that we cannot arbitrarily 
lower the quark mass at fixed lattice spacing, at least if we insist in integrating over 
continuum-like zero modes. 

Expanding the LO Lagrangian we obtain for the action the following decomposition 




(2.14) 




(2.15) 
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where 



S, 



(0) 



E 
2 



d xTr 



M n U + C^-M' 



[s< o) ; 






b 4 7 



d 4 xTr [^7r(a;)^7r(x)] 



5. 



(0) 



-ly'a 2 / (fx 



c 



Tr (c/ + Ut 



(2.16) 
(2.17) 



If we now use standard techniques we can write the partition function up to a normaliza- 
tion factor as 



Z = A d6 sin' 6e 



2 Z) ..21 COS 6 + 22 COS 



A / d6sm 2 6e zlCOse Y^(l 
Jo ^ n\ 



sin 



,2 mn 



where assuming a mass matrix proportional to the identity we define 

z x = 2m'EV, z 2 = 16a 2 VW. 



(2.18) 



(2.19) 



We see here the towers of infinite 0(a 2 ) terms induced by the coupling of the zero modes 
with the cutoff effects. Using the binomial formula we can write 



oo n n 

Zr, 



n= 

where we have defined the function 



n=0 U - k=0 



. k , 



-l) k X k+1 {zx) 



(2.20) 



VHXfc + 3/2) , , 



(2.21) 



and Ik+x{zi) are modified Bessel functions. To understand where this kind of expansion 
breaks down we perform the infinite volume limit at fixed lattice spacing and fixed quark 
mass trying to reproduce the same pattern of the results obtained in infinite volume [10]. 
So we want to study the behaviour of the partition function for large z\ and z 2 keeping 
their ratio fixed. The asymptotic expansion of the Bessel functions at leading order in 
1/zx reads 



so for large zx we have 

^£3E 



00 n n 



n=0 n - fc=0 



n 

. k , 



-1 



\ / 2irz 1 

fc 0Fr(fc + 3/2) e* 
{zx/2) k + l s/Wz~x 



(2.22) 



(2.23) 



The series can be rearranged in the following form 

2e zi+z 2 ~ f-2z 2 \ k F(k + 3/2) 



Z~A- 



2z\ k=o v 2 i 



fc! 



(2.24) 
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The radius of convergence of this power series with argument 



2 -: 2 



is given by 



k— >oo 



lim 



r(fc + 3/2) (fc + 1)! 
T(k + 5/2) k\ 



(2.25) 



which given the properties of the V function is 1. We conclude that if we take as a starting 
point the action in (2.15) at fixed lattice spacing and quark mass for large volumes the 
partition function is well defined only if 



We recall that the point Z\ = 2|z 2 | is exactly the location of the phase transition in the 
Aoki scenario. This implies that if we want to study the e regime with W%PT retaining the 
continuum-like integration over the zero modes we cannot lower the quark mass beyond the 
border of the phase transition. Moreover it is conceivable that given the strong coupling 
of the zero-modes with the cutoff effects it is better, to avoid large cutoff effects, to stay 
not so close to the phase transition point. 

The bound given in eq. (2.26) is obtained integrating over the continuum zero modes, 
thus it suggests that if we want to go below the bound we have to change the way we 
integrate over the zero modes. 

Equivalently in the infinite volume limit we are able to reproduce the well-known results 
concerning the chiral phase diagram of Wilson fermions [10] for quark masses not lower 
than the location of the Aoki phase transition. If we would like to go beyond this limit we 
would need to parametrize our integration over the zero modes taking into account the 
nature of the spontaneous symmetry breaking (SSB) taking place in infinite volume [21]. 

If we want to keep the description in terms of continuum- like zero modes one needs stay 
in the GSM regime or at most between the GSM and the LCE regime. 



3 Partition function and chiral condensate 

In this section we concentrate on the transition region between the GSM regime and the 
LCE regime. We try to understand the impact of the cutoff effects with Wilson fermions 
in the e regime if we simulate between the GSM regime and the LCE regime treating the 
operators parametrising the cutoff effects as insertions. 

From the power counting point of view this can be seen in two ways, either as performing 
the computation up to 0(e 2 ) using the GSM power counting and including terms of 
0(e 4 ) containing the relevant operator parametrising the cutoff effects, or performing the 
computation up to 0(e 2 ) in the LCE, but treating the 0(a 2 ) cutoff effects terms as small 



Zx > 2\z 2 \, 



with 



zx > 



(2.26) 
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compared to the mass term. In both cases the pattern of spontaneous symmetry breaking 
is as in the continuum. In our analysis we want to keep 0(a 2 ) and 0(am) terms. Keeping 
the 0(a 2 ) term is equivalent of retaining the first term of the expansion done in the 
previous section which interpolates between the two regimes. We add also the O(am) to 
understand the impact of this term in the correlation functions. 

We have already seen in sec. 2.1 that the LO Lagrangian in the GSM regime is like in the 
continuum 



d ll U{x)d li U( X y 



Tr 

2 



M n U{x)+M'U{ 



x 



(3.1) 



as the pattern of spontaneous symmetry breaking. We can thus use the standard parametri- 
sation for Goldstone fields (2.12). Expanding the LO Lagrangian we obtain for the action 
the following decomposition 



B 



(3.2) 



where 



S, 



(o) 



A,B 



are defined in eq. (2.16), which are the standard LO continuum contribu- 



tion from the zero modes. 

At the next order we have two types of contributions coming from the expansion to 0(e 2 ) 
of the LO Lagrangian and from the cutoff effects. We write then the next term of the 
action as 



s 4 = si 2) + ss. 



where 



B 



is the standard Ofe 2 ) contribution in the continuum 



S. 



(2) 



a 48F 2 
E 



d 4 xTr {[8^71, 7r] [O^tt, it]} 



b AF 2 



d xTi 



M'Uqtt 2 + iPulM' 



while 



(3.3) 
(3.4) 

(3.5) 
(3.6) 



'„2 



5S A = J d\ 
2aWE 



Tr (U + U ] )f + ^Tr [d^d^] Tr [u + U ] Q 



F 2 



-Tr 



M'm + W M' Tr U + U { 



F 2 



Tr 



M n + M! 



(3.7) 



These are the terms we want to include containing only lattice artifacts. The last term 
is responsible for the contact term in the chiral condensate. In the following we will call 
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the inclusion of and SS^ an 0(e 2 ) expansion keeping in mind that includes few 0(e 4 ) 
terms in the GSM power counting. 

The measure of the partition function can be splitted with appropriate Jacobian J (tt) 

Z = Jv[U] e~ s[u] = J V [U ] V [tt] J (vr) e~ s[u] (3.8) 
and written in the following exact form 



Z= V[U ]e 



AO) 



^ [U 



where 



S[U]- 



? (o) 



(3.9) 



(3.10) 



Z V [U ] = jV[n]J(n)e 

To expand the partition function up to 0(e 2 ) we need to expand the Jacobian and the 
exponential in eq. (3.10) up to 0(e 2 ). For convenience we rewrite the LO partition function 

as 



= Z®Z<?\ where zf ] = J V [U ] e 



; (0) 



2 (0) 



j(0) 



4 2) 



B 



(3.11) 
and terms 



Expanding to 0(e 2 ) [Uq] we obtain terms which depends on Uq like 
that, like the Jacobian, do not. These latter will enter only in the absolute normalization 
of the partition functions. We then end up with several terms 



HV (n 2 - 1 

zAUo]=^{i U^ Tl 



M^U + UlM] G(0)+ 



+ W'a 2 V 
2aZH' 



+ 



F 2 



-Tr 



Tr (U + Ul 
M' ] + M' 



2 laWYM 



Tr 



M^Uq + U$M] Tr \U + Ul 



+ 



(3.12) 



where G(x — y) is the 'pion' propagator without the contribution of the zero mode defined 
by 



7 (o) 



V [tt] e" 



j(0) 



'7r„(z)7r fc i(y) = 2 U> a 5 kl - -^-^ijhij G(x - y). 



In dimensional regularization the propagator G(0) is finite and is given by 



G(0) 



with V = L 6 x T 



(3.13) 



(3.14) 



and (3\ is a numerical constant which depends only on the geometry of the box [11]. 
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In eq. (3.12) the first term that can be reabsorbed in the definitions of S, thus defining 
an effective LEC £ e fr = Sp where 



p = l 



N2 f 



N f F 2 



1 G(0). 



(3.15) 



The partition function at 0(e 2 ) can now be written as 



Z=M I V[U ]e 

2aWEV rT1 
-Tr 



,(0) 



Peff] 

a x i 1 + W a z V 



F 2 



M j U + uIm] Tr \U + U { 



+ 



2aLH' . 



F 2 



-Tr 



+ 



M' ] + A-f' 



(3.16) 



We can now analyze various features of this formula. The first thing to notice is that 
with the choice of our power counting the Boltzmann factor is as in the continuum, while 
cutoff effects appear as standard perturbative corrections. In the continuum the 0(e 2 ) 
contribution to the partition function is reabsorbed in the definition of E. Here at finite 
lattice spacing we define have additional terms that cannot be reabsorbed. 

We can now proceed and compute the chiral condensate. It is well known that the direct 
computation of the chiral condensate is a difficult task with Wilson fermions. Here we 
simply use it as an example to explain the computational procedure. We first define the 
quark currents in a standard fashion 



S°(x) = ip(x)i/j(x) 



P^x) 



i){x)i^tl)[x) 



(3.17) 



S a (x)=^(x)-^(x), 
where ip is a flavour doublet quark field. 



T 



P a (x) = i)(x)i^—i)(x) 



(3.18) 



To compute correlation functions in the effective theory we use the standard method of 
augmenting the mass term with sources for scalar and pseudoscalar currents 



M 



M. + s(x) + ip(x) 



(3.19) 



where 



six) 



s°(x)+s a (x 



T 



T 

p(x) —p(x) +p a (x) — . 



(3.20) 



The corresponding correlation functions can now be obtained performing functional deriva- 
tives respect to the partition function (3.9), which now depends on the sources, and then 
setting the sources to zero. 
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We can now write the scalar and pseudoscalar current at NLO in the effective theory 



- S°(x) 



-Tr 



U(x) + U(x) ] 



Tr (U(x) + U(x) n 



F 2 



F 2 



(3.21) 



P a (x) = -Tr 
v ; 2 



IT 



— {U(x) - U{xf 



2a£W m 


-%T a , 


+ Tr 




P2 





Tr (U{x) + U{x) ] 
(3.22) 

Expanding this expression to 0(e 2 ) and performing the contractions on the tt field, in the 
same way as we did for the partition function, we obtain for the chiral condensate the 
following expression 



(S°) 



eff (Tr U + Ul) 2 + 



([Tt(u + U^] 2 )o + 



where we have defined two different expectation values over the zero modes 

1 



(3.23) 



(O[U ]) 



z, 



(0) 



V[U ]O[U ]e 



AO) 



(3.24) 



and 



(O [U }) 2 = (O [U ])o - (O [U ] 5S A [U })o + (O [U ]) (8S 4 [U ]) . (3.25) 

In the last two terms of this formula we can compute the expectation values not consid- 
ering the effective LEC S e fj but only E, the difference being of higher order. We define 
now the following LECs 

lew , w 

6 w = ~pr, w = -yr- ( 3 - 26 ) 

This definition is similar to the one done in ref. [17], but here we keep explicit the a 
dependence in the correlation functions and we redefine the LECs absorbing the Wq term. 
As a result their dimension is different and given by [5w] = [Energy] and [it/] = Energy 4 . 
These LECs have a very precise origin: 5w parametrises the mass dependent 0(a) cutoff 
effects of the pion mass and it/ the 0(a 2 ) cutoff effects of the pion mass. Setting c^w to 
its 'correct' non-perturbative value sets 8w = but simply changes the value of w* 2 



The final expression for the chiral condensate reads 




(3.27) 



2 If we use Wilson twisted mass fermions the same LECs play a relevant role: 5w parametrises 
the dependence of the PCAC mass on the twisted mass close to maximal twist and w' the 0(a 2 ) 
pion mass splitting. 
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Where z 2 = a 2 w'F 2 V, the argument of the functions is z = 2mS c fjV A and we have defined 



7T i7T 

X 1 (z) = -h(z), X 2 (z) = —I 2 (z) (3.28) 



z z 2 



where I\ and I 2 are modified Bessel functions (definitions and properties can be found 
for example in [22]). Important checks are obtained performing the infinite volume limit 
and the continuum limit. If we send o -> we obtain the 0(e 2 ) result of Gasser and 
Leutwyler [1]. If we send V — > oo we obtain a formula for infinite volume and finite 
lattice spacing. This is consistent with the formula given for example in ref. [17], if we 
just compare the terms we have actually included in our computation. In particular the 
terms which are proportional to a 2 V or az both vanish in the infinite volume limit because 
for large z 

In the left plot of fig. 1 we show the dependence of the chiral condensate on the quark 
mass, for a lattice spacing given by z 2 = 1 and a volume of V = L 3 x 2L with L = 2 
fm, where we have subtracted the contact term. We also use the values for the LECs: 
£ = (250MeV) 3 , F = 86.2Mev, ad w = 0.1. 



4 Two-point functions 



In the previous two sections we have explained in detail the strategy we wish to adopt, 
and we can repeat a similar computation for scalar and pseudoscalar two-point functions. 
We define 

-i J d 3 x(S°(x)S°(0)) = C s (z ), ^ / d 3 x(P a (x)P b (0)} = 5 ab C P (x ), (4.1) 
and after performing the group integrals we obtain (X 3 (z) is defined in eq. (2.21)) 



C s (x )=4E e 2 fr 1 



Xo(z) 



XAz) 



+ (z 2 + aSwz) 



XJz) 



X 2 (z) 



Xi(z) X 2 (z) 



+ 2a5w 



X[(z) 



X' 2 {z) 



X 1 (z) 



+ 



8aH'X[( 

T^X^iz) ' F 2 X l (z)L 3 ' 



(4.2) 



where r = Xq/T and where the time dependence of the massless propagator is given, as 
in the continuum by [11] 



hi(r) 



1 

f 



d 3 xG{x) 



T 
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(4.3) 



For the pseudoscalar correlation function we obtain 
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z 2 = 0.5 NLO 



m D [MeV] 



15 



3C p (x )/E eff 




0.2 0.4 0.6 



Fig. 1. Left plot: mass dependence of the subtracted and properly normalized chiral condensate in 
the continuum and at finite lattice spacing. The dashed lines indicate the range of quark masses 
where the expansion carried out in this work breaks down. This range shrinks to zero towards the 
continuum limit. Right plot: Euclidean time dependence of the pseudoscalar two-point function 
in the continuum and at finite lattice spacing with 2 values of a5\y = 0,0.1. 



C P (x ) 




z 2 + a5y/z) 



Xi(z) 



X 1 (z) X?(z)_ 



3 / lX 2 (z)\ T 



(4.4) 



Using the same values for the LECs' we have used for the chiral condensate in the right 
plot of fig. 1 we show the Euclidean time dependence of the two-point functions. 



Formulas (4.2) and (4.4) suggests an interesting way to remove most of the cutoff effects, 

which is to study the correlation function Cs(xo)/4 + 3Cp(xo). In fact in this correlation 

r x' (z) i 

function all the cutoff effects are absent but the term proportional to 2a5w x\(z) ' ^ e 

can immediately recognize the ratio xT^j ^ as ^ e con tinuum mass dependence of the 
chiral condensate in the symmetry restoration region. So this term is expected to be very 
small in the e regime and can be removed determining in a non-perturbative way csw- 
The correlation function will still have a suitable dependence on the LECs, so it can 
still be used to fit lattice data. The obvious drawback is that Cs(xo) would require the 
computation of a disconnected diagram. 
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5 Conclusions and outlooks 



In this Letter we have analyzed some aspects of the chiral limit with Wilson fermions. We 
have discussed how cutoff effects influence the way we reach the e regime, and to which 
extent we can continue to discuss about zero-modes as it is done in the continuum. We 
have found that concerning cutoff effects there is an abrupt transition between the GSM 
regime, where Wilson fermions show no cutoff effects up to NNLO, and the LCE regime 
where cutoff effects appear at LO directly in the action, in a non-perturbative way. This 
transition is driven by the coupling of the zero modes with 0(a 2 ) effects. A first natural 
way to avoid these cutoff effects is to keep the quark mass heavier than the values at 
which one would expect phase transitions in infinite volume, or better in the GSM regime 
where Wilson fermions are automatically 0(a) improved up to NNLO and continuum 
formulae could be used. One can also simulate in the transition region between the GSM 
and the LCE regimes. For those values of quark masses we can still use W^PT in the e 
regime in a very similar way to what is done in the continuum, i.e. keeping the integration 
over the SU(2) manifold for the zero modes, and treating the cutoff effects as perturbative 
corrections. In this transition region we have computed formula? which describe the impact 
of the cutoff effects as soon as we move away from the GSM regime. This formulae could 
be used to fit lattice data in order to better understand the systematics of the simulation 
results coming from discretization errors. 

The analysis presented in this Letter indicates a second way on how to control the increase 
of the discretization errors in the transition from the GSM and the LCE regime. Our 
results give strong indications that the rise of the cutoff effects is driven by the zero 
modes. It is possible that, if one is able to remove the zero modes from the analysis of 
the lattice data, the cutoff effects will be again under control. This could be done for 
example finding a way to define in a robust way configurations with trivial topology. The 
analysis presented here seems to indicate that the cutoff effects blow up only in sectors 
with non-vanishing topological charge, i.e. with zero modes. 

Note added 

After this work was completed and sent to the whole ETMC I was informed and then 
I received a paper [23] about Wilson fermions in the e regime where similar conclusions 
have been reached. 
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